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ABSTRACT 

Context. To date the onset of large-scale star formation in galaxies and its link to gravitational stability of the galactic disk 
have not been fully understood. The nearby face-on spiral galaxy M51 is an ideal target for studying this subject. 
Aims. This paper combines CO, dust. Hi, and stellar maps of M51 and its companion galaxy to study the H2/H1 transition, 
the gas-to-dust ratios, and the stability of the disk against gravitational collapse. 

Methods. We combine maps of the molecular gas using ^^CO 2-1 map HERA/IRAM-30m data and Hi VLA data to study the 
total gas surface density and the phase transition of atomic to molecular gas. The total gas surface density is compared to the 
dust surface density from 850 fim. SCUBA data. Taking into account the velocity dispersions of the molecular and atomic gas, 
and the stellar surface densities derived from the 2MASS if-band survey, we derive the total Toomre Q parameter of the disk. 
Results. The gas surface density Egas in the spiral arms is ~ 2 — 3 higher compared to that of the interarm regions. The ratio 
of molecular to atomic surface density shows a nearly power-law dependence on the hydrostatic pressure Phydro- The Egas 
distribution in M51 shows an underlying exponential distribution with a scale length of hgas ~ 7.6 kpc representing 55% of the 
total gas mass, comparable to the properties of the exponential dust disk. In contrast to the velocity widths observed in Hi, 
the CO velocity dispersion shows enhanced line widths in the spiral arms compared to the interarm regions. The contribution 
of the stellar component in the Toomre Q-parameter analysis is significant and lowers the combined Q-parameter Qtot by up 
to 70% towards the threshold for gravitational instability. The value of Qtot varies from 1.5-3 in radial averages. A map of Qtot 
shows values around 1 on the spiral arms indicating self-regulation at play. 

Key words. Keywords should be given 



1. Introduction 

To understand the evolution and appearance of galaxies 
it is crucial to study the interplay between stars and the 
interstellar gas and dust. Stars are made within the dense 
and cold interiors of molecular clouds. It is generally be- 
lieved that molecular clouds form from the densest regions 
of more widely distributed, diffuse atomic clouds. During 
their lifetimes, stars seed the gas with heavy elements, en- 
ergy, and momentum, all of which strongly affects their 
parental environment, the subsequent formation of new 
stars, and their host galaxy. The photo-destruction of 
molecular clouds leads to their dispersal and the forma- 
tion of diffuse atomic clouds, thereby recycling the ma- 
terial. The processes governing the phase transition from 



the cold neutral medium to dense, molecular gas, and back 
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again, are however still unclear (IRosolowskv et al.l 12007 



Knapen et al. 20061 : Heitsch et all 2008 ). Studies combin- 
ing maps of CO, the traditional tracer of molecular gas, 
Hi, and the cold dust will form the basis for any under- 
standing of their interrelation. The spatial distribution of 
the various galactic components is best studied in nearby 
face-on galaxies, like M51, w hich lies at a distance of 



4Mpc (jFeldmeier et al.lll997t) 



Two empirical laws describe star formation in galaxies. 
First, stars are observed to only form efficiently above a 
critical gas surfa ce density tha t appears to be determined 
by the Toomre ( Toomr j 1964 ) criterion for gravitational 
stability. This star for mation threshold has been o bserved 
in many galaxies (e.g. Martin fc Kennicutd 200l | ). Note- 
howev er the exceptions discussed e.g. by IWong fc Blitd 
(|2002l ). Second, a tight relation has been found between 
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Fig. 1. Left: Map of integrated CO 2-1 intensities on the antenna temperature scale in Kkms ^ at 11" resolution 
(cf. Paper I) . A loga rithmic spiral is shown in all Figures of M51 to guide the eye. Right: Map of Hi column densities 
([Walter et al.ll2008h in units of 10^° cm'^ at a resolution of ~ 6". 



the star formation rate and the surface density of the total 
gas, i.e. the molecular and atomic gas, the Schmidt law. 
In this paper, we will focus on the Toomre criterion. 

In its original version it considers a dynamical stability 
analysis of a thin, gaseous, differentially rotating, single- 
component disk against axisymmctric gravitational per- 
turbations. Observational studies of the Toomre threshold 
have usua Uy assumed a constant velocity dispersion of the 
gas, while Schavd ( 2004[ ) questions this assumption, argu- 
ing that the drop of velocity dispersion associated with 
the transition from the warm to a cold phase of the inter- 
stellar medium causes the disk to become gravitationally 
unstable. In addition, recent studies have shown that the 
stellar component contributes significantly to the gravi- 
tational stability and c annot be neglected ( Boissier et al.l 
2003HYang et al.ll2007t ). 



In I Schuster et al.l (|2007l .hereafter Paper I), we com- 
bined new CO 2-1 data of M51 obtained at the IRAM- 
30m telescope using HERA, with H i, and radio continuum 
data from the literature to discuss radial profiles of the 
molecular and atomic gas mass surface densities and star 
formation rates, averaged over azimuth as a function of 
radius in the disk of M51. In the following, we will simply 
speak of radial averages. 

In the present paper, we combine the CO 2-1 
data with Hi VLA data from the THINGS survey 



(IWalter et allbOOSl l2008l) SSOu m SCUBA dust contin- 
uum data (jMeiierink et al.ll2005l ). and a if-Band 2MASS 
image. We use the CO and H i data to derive maps of the 
mass surface densities of the molecular and atomic gas, 
and maps of their velocity dispersions. Maps of the total 
gas and of the ratio of molecular to atomic gas are then 
discussed. We separate the total gas mass surface density 
map into an underlying exponential disk and the spiral 
arm structure. The results are then compared with the 
dusty exponential disk found bv lMeijerink et all I 2005[ ) in 
their SCUBA map of M51. This allows us to study the 
variation of the gas-to-dust ratios in the disk and in the 
arms. 

In preparation for the Toomre analysis, we discuss 
the velocity dispersions of the molecular and atomic gas. 
Combining the data of the molecular and atomic gas with 
the map of stellar mass surface density derived from the 
2MASS data, we then present the radial average and a 
map of the total Toomre Q parameter describing the grav- 
itational stability of the disk. 



2. Observations and Data 

Observations of the ^^CO 2-1 emission from M51 were 
conducted with the IRAM 30m telescope in February 
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Fig. 2. Map of the total gas mass surface density Sgas 
in M51 in Mq pc"^ derived from the ^^CO 2-1 and Hi- 
emission. The cross marks the (0,0) center position at a = 
13^'29"52!7, S = 47°11'43" (eq.2000). The 1-cr value of 
the total gas density map is 0.8 M0 pc~^. Contours mark 
regions discussed in ij3.ll and Tabic [1] 



2005. The half power beam width (HPBW) is 11". For de- 
tailed information of the observations we refer to Paper I. 

Additionally, w e use data of the Hi N earby Galaxy 
Survey (THINGS) (|Walter et al.ll2005ll2008f ). M51 was ob- 
served in the VLA D, C, and B arra y configurations and 
presented in iKennicutt et al ( 2007 ). In D-configuration 
the interferometer is sensitive to scales of 15' allowing to 
detect almost the total flux for M51. The angular reso- 
lution in the integrated intensity image is 5.86" x 5.56". 
The sensitivity of the THIN GS data is 0.44 mJ y beam^i 
for 5.2 km/s channel width ( Walter et al. 2008|). which is 
about a factor 4 deeper than the VLA data of [Rots et aT 
( 1990f) used in Paper I. 



To derive stellar densities, we used the K-Band images 
of the Two Micron All Sky Survey (2MASS) Large Galaxy 
Atlas (jJarrett et al.ll2003h . 



3. Molecular and atomic gas 

3.1. The total gas surface density 

In Fig[2]wc present a map of the total gas surface density 
Sgas. This map is a major ingredient of the gravitational 
stability analysis we will discuss later in this paper. 



To convert from the observed integrated intensities 
to the total gas surface densities, we use the relation 
5^gas = 1-36 X (Sh2 + SIhi) taking into account the mass 
contribution from He (cf. Paper I). The atomic surface 
density was calculated from the THINGS H i column den- 
sity image of M51 (Fig[T]), smoothed to 11" resolution. 
Two logarithmic spirals (shown in black) serve as a guide 
to delineate the inn er spiral arms. They are adopted from 
Shettv et all (|2007t ) who identified them from a BIMA 
i^CO 1-0 map of M51. 

The total gas surface density in Figl2] peaks at 127 
M0 pc~^. The inner spirals showing the densest peaks just 
south-west of the center show total gas surface densities 
of 80-100 Mq pc-2. 

We subdivided the total gas map (Fig.[21) into five 
different regions: the inner and outer spiral arms, the 
intcrarm regions, the outer disk, and the companion 
galaxy. All regions above a threshold surface density of 
20 Mq pc~^ were attributed to the spiral arms. The outer 
border of the outer disk of M51 was defined by the 3ct 
level of (Tgas = 2.4Mopc-2 ((t(S]h2)= 0.36 Kkms-\ 
(t(Shi )= 0.44 Kkms~^). The threshold for the arms 
leads to the identification of spiral arms which coiie- 
spond roughly with the a rms seen in interferometric CO 



maps IShettv et al.l (|2007r ). and thos e seen in Ha-, and 



in 20 cm radio-contin uum emission (jScoville et al, 
Patrikeev et"ani2006f) . 



2001 



Table[T] lists the average total gas density and the av- 
erage H2 over Hi density ratio for the five subregions. 
The average of the inner spiral arm region, is Sgas ~26.8 

Mq pC-2. 

The outer spiral arms show less gas surface density, 
with an average total gas density of 18.8 Mq pc~^, includ- 
ing the south western spiral arm and the north eastern 
spiral arm. This is approximately a factor of 1.5 weaker 
than the spiral arms in the central region. Few high den- 
sity regions on the outer spiral arms reach Egas of 40- 
50Mqpc-2 e.g. to the south-west at (-120", -200") and 
the north-east at (140", 100"). This is signifcantly weaker 
compared to the peaks in the center. 



Table 1. Properties of Egas and Eh2/Shi in M51 averaged 
over the areas shown in Fig[2] 



[M© pc 



gas 
-21 



Sh2/Ep 



inner spiral arms 


26.8 


2.3 


outer spiral arms 


18.8 


0.21 


interarm region 


9.4 


0.71 


outer disk 


6.4 




NGC 5195 - companion 


20.2 


0.4 


ratio inner/outer spiral 


1.4 


11.0 


inner arm/interarm 


2.9 


3.2 


outer arm/interarm 


2.0 


0.3 
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radial average (cf. Paper I). The spiral arms in the center 
are the regions with the highest fraction of molecular gas. 
The average ratio in the two inner spiral arms delineated 
in Fig[2]is 2.1 (Tablc[T|). Outwards, at radii of more than 
about 100", the ratio drops also in the spiral arms. The 
north-eastern and south-western continuation of the inner 
spiral arms show few regions with a ratio of 1 or slightly 
above which are correlated with the regions of increased 
total gas density (Fig[21). The average ratio on the outer 
spiral arms in the marked polygons is 0.21. This is a factor 
of 10 lower compared to the inner spirals. The interarm 
region within the inner 100" radius shows an average ratio 
of 0.71. The molecular gas surface density drops below the 
2(T threshold at larger radii for most positions. 

The increase of the molecular-to-atomic gas fraction 
in the spiral arms may be triggered by the spiral wave 
increasing the hydrostatic pressure, and thereby the for- 
mation of molecular clouds from the atomic phases of the 
ISM. 

Next, we will discuss the hydrostatic pressure in the 
disk of M51 as a possible parameter governing the ratio 
of molecular and atomic gas. 

3.3. The hydrostatic pressure 



Fig. 3. Map of the ratio of molecular- (SH2) a-nd atomic 
gas surface density (Shi) m M51. Pixels below 2a of the 
masked- moment integrated intensity image of ^^CO 2- 
1 are blanked. The l-cr limit of the ^^CO 2-1 data is 
0.65Kkms~^ corresponding to Sh2= 0-36 Mqpc"^. As 
the Hi is more extended, the threshold of the molecular 
gas limits the extend of the map. 

The interarm regions show surface densities of 5-15 
M0 pc~^. The average value is 9.4 M0 pc~^. The averaged 
inner arm/intcrarm contrast in the t otal gas surface den- 



sity is 2.9 confirming the results of iGarcia-Burillo et al.l Pj. 



(19933) in the integrated CO 2-1 intensities. 

There are several spur-like structures bridging the gap 
between the inner spiral structure and the outer spiral 
arms e.g. in the north-western area at (—40", 80") at a 
surface density of 20-30 M0pc~^. Also, in the north-cast 
at (100", 100") there is a strong bridge between inner and 
outer spiral arms with a total gas surface density of 20-40 
Mq pc~^. The gaps or minima in total gas density seen in 
the spiral arms e.g. at (-80", -40") and at (-120", 60") 
are interpr eted as the signature of 4:1 ultra-harmonic res- 
onances bv lElmegreen et al ( 1989[ ) applying density wave 
theory. 

3.2. The molecular gas fraction 

Figure[3] shows the ratio map of molecular versus atomic 
gas surface density Sh2/Shi. Overall, the map shows a 
steep decrease of the ratio from peak values of around 50 in 
the center to 0.1 in the outer regions as already seen in the 



Recent invest igations, e.g. by iBlitz fc Rosolowskv 
(|2004al |bl. I2OO6I ). have shown a tight correlation of the 
Eh2/Shi ratio i n external g alaxie s and the hydrostatic 
pressure Phydro- lElmegreenl (1991) theoretically studied 
the impact of hydrostatic pressure and radiation field on 
the ratio of atomic and molecular gas. He concludes that 
the pressure should be the dominant factor determining 
this ratio. He predicts a powe r law ~ f of thei r 
radial averages. Following iBlitz fc Rosolowsicvl (HooiB), 
the hydrostatic pressure can be estimated for an infinite 
disk with isothermal gas and stellar layers from the 
midplanc pressure in equilibrium: 



hydro 



(2G) 



0.5 



^gas'^gas 



(1) 



Assuming the volume density of the gas pg^s is small com- 
pared to the stellar density and a self-gravitating stellar 
disk (E, = 2p,/i»): 



hydro 



0.84(GI],)°-^I]g 



'^gas 



(2) 



I]» is the stellar surface density, which we will address 
in a later section, Sgas the total gas density and o-gas the 
velocity dispersion of the gas, which we estimated from the 
^^CO 2-1 velocity dispersion (cf. Paper I). For the stellar 
scale height, we adopt hn, = Ikpc, c onsistent with stud - 
ies bv lvan der Kruit fc Searld (|l98ll ): lKregel et all (|2002l ) . 



-Phydro in equation [5] is given in Nm~^, assuming all input 
parameters are in Sl-units. 

As in Paper I, radial averages were created by averag- 
ing in elliptical annuli spaced by 10", including and taking 
into account points of no detection, and centered at the 
center position. The annuli arc circular rings viewed at an 
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Fig. 4. The radial averaged fraction of molecular gas den- 
sity and atomic gas density Shi is plotted against 
the hydrostatic pr e ssure Phydro- The slope found in 
Blitz fc Rosolowskvl (|2006l) is indicated in red In green 
the linear regression fit to our data P^* is shown. Going 
from left to right galactocentric radius decreases. 



inclination angle i = 20° and with a line of nodes rotated 
from north to east by the position angle PA=170°. 

FigureH] shows the SHa/^Hi ratio plotted versus the 
hydrostatic pressure Phydro, covering three orders of mag- 
nitude in pressure, the largest range measured so far in a 
galaxy disk. A linear regression fit yields a powerlaw: 



^H2 



(fliydro/Po)", 



(3) 



with the powerlaw coefficient a — 0.87 ± 0.03 and Pq= 
(1.92 ± 0.07) 10^) cm-^K. The slope of 0.87 is very simi- 
lar to the slope of 0.92 determined bv lBlitz fc Rosolowskv 
( 2006h for a sample of 14 spiral and dwarf galaxies, in- 
cluding M51. We find a slightly lower Pp compared to the 
mean result of iBlitz fc Rosolowskvl (|2006f ) in their sample 
with (3.5 ± 0.6) 10^) cm'^K. 

This result underlines the importance of the hydro- 
static pressure as the physical parameter determining the 
fraction of molecular to atomic material at a given ra- 
dius on large scales. The nearly linear dependence found 
over the wh ole pressure regime is expected by theoreti cal 
predictions (|Blitz fc Rosolowskvl2006HElmegree^ll993l ). if 
the gravitational potential imposed by the gaseous com- 
ponent is small compared to the stellar component. Note 
that the pressure dependence of the molecular to atomic 
fraction is expected to break down when the spatial res- 
olution of the observational data is close to the size of a 
single giant molecular cloud (GMC). A single CMC is ex- 
pected to have enhanced pressure due to self-gravity and 
the assumption of pressur e equilibrium is no longer valid 
( Blitz fc Rosolowskvlliooeh . 
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Fig. 5. The residual map of the total gas surface density 
after removal of a fitted underlying exponential disk. The 
exponential disk is shown in contours with contour levels 
IMq pc~^, 5M0 pc^^ increasing to ISM© pc~^ in steps of 
5M0 pc-2. 



3.4. An exponential disk in M51 



Table 2. Properties of the exponential disk of the total gas 
derived in this paper compared to the properties of the expo- 
nential dust disk determined bv lMeiierink et al] |200^. 



Total gas exponential disk: 




Egas scale length h^^ 


7.6 ±0.1 kpc 


Amplitude: E^if 


17.8 ±0.2 Mopc-^ 


disk/total mass-fraction gas 


61% 


Dust exponential disk: 




dust scale length hd 


6.3 kpc 


Amplitude: T.^^'J^ 


0.77 ±0.2 Mopc-^ 


disk/total mass-fraction dust 


55% 



Next, we investigate the model of an underlying expo- 
nential distribution of the total gas density in the disk 
of M51 and compare it with the dust disk found by 
Meiierink et al. lilooi). We fitted an inclined elliptical ex- 
ponential disk to the Egas total gas surface density dis- 
tribution map. We assumed a centered disk at an inch- 
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Fig. 6. The residual map of the 850/ini dust emission after 
removal of a fitted u nderlying exponential disk (Fig. 2 in 
Meijerink et al.ll2005h . The Sa-level is at 24 Jy beam"^ 



nation of 20° and a position angle of 170°. The scale 
length hgas and the amplitude S|as'^ are determined from 
a 2-d fit of an exponential distribution Egas ~ T,p^^^x 
exp(— i?//igas). AH pixels in the total gas density map ex- 
ceeding a threshold surface density of Egas = 20 M© pc~^ 
were blanked. This corresponds to the arm regions we 
marked in Figl2] and discussed in A significantly 

higher blanking threshold leads to an overestimation of 
the exponential disk in the central parts as spiral arm re- 
gions are included in the fitting. 



Properties of the exponential disk The scale length of the 
fitted exponential decay of the total gas disk is /igas — 
7.6 kpc (Tabled!) with the peak at 17.8 Mepc'^. The 
reduced of the fit is 5.1. The fraction of total gas mass 
contained in the exponential disk is 61%. Thus the spiral 
ar ms contain only ^ 40% of the total gas mass. 

Meijerink et al. ( 2005[ ) fit a scale-length of the ex- 
ponential dust disk of /id = 5.45 kpc corresponding to 
hd = 6.3 kpc for a distance of 8.4 Mpc, assumed in the 



present paper. Assuming an isothermal disk, the fraction 
of dust contained in the exponential disk is 55%. Notably, 
the scale lengths of the total gas and of the dust exponen- 
tial disks are similar, as are the relative fractions of mass 
contained in the disks. 

Note that the shallow radial temperature profile de- 
duced by iMeiierink etldl (|2005l ) of 25 K in the center to 
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Fig. 7. (a) Radial averages of surface densities. E* is the 
stellar surface density. Drawn lines show averages includ- 
ing the companion galaxy. The dashed line shows stellar 
densities, ignoring the companion. Egas is the surface den- 
sity of the total gas. E™s is the density of the fitted ex- 
ponential gas disk. Edust is the density of the dusty disk, 
^dust ^^'^ density of the exponential dusty disk. Dust 
values below 3a = 24 Jy beam~^ are ignored, (b) The gas- 
to-dust ratio for the exponential disk and for the residual 
data, i.e. the spiral structure. 



17 K in the outskirts at 10 kpc radial distance, only leads 
to a 35% variation of corresponding dust surface densi- 
ties if the temperature is assumed to be constant. Thus, 
an isothermal dust disk of 25 K dust t emperature and 
'^850 = 1.2cm^g~^ (jMeiierink et al.ll2005l ) are assumed in 
the calculation of the dust surface density. The dust sur- 
face density Edust is 



-'dust 



c 



i^850dV('«850B850(Tdust)). (4) 



-Fsso = 'S'85o/5^850 X ^beam is the integrated flux. S'gso is 
the flux density, fisso the beam size of 15" and Abeam the 
beam area at the distance c?=8.4 Mpc. Ssso labels the 
Planck- Function at 850 fim. 

The radially averaged gas-to-dust mass ratio of the 
surface densities of the exponential disk is shown in the 
lower box of Figurc[71 The gas-to-dust mass ratio is nearly 
constant with galacto-ccntric radius, reflecting the similar 
scale lengths. Values vary only between 23 and 26, which is 
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about a quarter of the canonical Galactic gas-to-dust mass 
ratio of 100. This constancy is in contrast for example to 
the strong radial variation of the H2/H1 surface density 
ratio by more than a factor of 100. 

The gas-to-dust mass ratio in the exponential disk 
is significantly lower compared to the canonical Galactic 
value, whereas the metallicity of M51 is only slightly su- 
persolar with a s hallow gradient with varying galacto- 
centric radius (e.g. iKramer et al.ll2005l ). However, several 
sources of uncertainty enter into this calculation which can 
not be resolved in the context of this study but might drive 
the value towards the canonical gas-to-dust mass ratio of 
100. In the inner part H2 dominates the total gas mass 
and is itself calculated via the X-factor for M51 (see dis- 
cussion in Paper I) . The low gas-to-dust mass ratio in the 
exponential disk, which constitutes the interarm regions, 
might indicate an X-factor in the interarm regions closer 
to the Galactic X-factor (4 times the value of M51). This 
will on the other hand not severely effect the gas-to-dust 
mass ratio at larger radii where the atomic component is 
dominanting the molecular contribution. 

Additionally, the assumption of a constant dust 
temperature Td i.st = 25K and Ksnn = 1.2cm^g^^ 



(jMeiierink et al.ll2005f) might cause an overestimation of 
Edust in the exponential disk. A 4 times higher Kgso yields 
a Galactic gas-to-dust mass ratio value but seems un- 
likely as the values for k^^q determined in ex ternal galax- 
ies a r e typically khro = 0.7 — 2.1cm^ g~-^ (I James et al. 



l2002t I Alton etaL 2002[ ). A higher dust temperature also 
drives the gas-to-dust mass ratio towards the Galactic 
value, with a ratio of 100 for Tjust = 121K and Kgso = 
1.2cm^g~^ as used above. 

3.5. Properties of the residual emission 

Figures[5]and[n]show the residual total gas and dust maps, 
after subtracting the corresponding exponential disks. The 
lower box in Figurc[7] shows the corresponding radial av- 
erages of the gas-to-dust mass ratio. 

For the computation of dust surface density in the 
residual data, i.e. in the spiral arms, we assume a uniform 
dust temperature Tdu st = 25K and ksso = 1.2cm^g^^ 
( Meiierink et al. I I2OO5I) as above. 

The resulting fraction of Egas and EJj^"j. in the spiral 
arms, shown in Figure[71 is not constant with radius, but 
rather varies by almost a factor 3 between about 70 and 
190 within the first 5kpc where the emission is above the 
3(T-threshold. 

In contrast to the gas-to-dust fraction in the exponen- 
tial disk, we find a minimum of 73 at 3.8 kpc radial distand 
and a value of 186 in the center. Thus the gas-to-dust mass 
surface density ratio is a factor of 3-8 higher in the spiral 
arms compared to the underlying exponential disk. 

As the exponential disk constitutes a minor fraction 
of the gas mass in the spiral arms, it hardly mixes with 
the spiral arm component. 



The global gas-to-dust mass ratio, summing the mass 
components in the exponential disk and the spiral arms in 
M51 upto the 3cr-limit of the dust observations, is at 112 
close the Galactic value. 



4. Stellar surface density 
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Fig. 8. Map of the if -band luminosity Lk from the Two 
Micron All Sky S urvey (2MASS) Large Galaxy Atlas 
( Jarrett et al.ll2003h smoothed to 11" resolution. The imits 
are Lq. 



The stellar surface density is an important parameter 
in determining the stability of the disk. We use the K- 
Band images of the Two Micron All Sky Survey (2MASS) 



Large Galaxy Atlas (jJarrett et al.ll2003f l to determine the 
stellar surface mass density E*. The if— Band is a reason- 
able tracer of the stellar mass density as it is much less 
affected by extinction within M51 than e.g. the i3— Band 
in the optical. The Galactic foreground extinction in the 
if— Band is 0.013 mag comp ared to 0.152 mag in the 



.B-Band (ISchlegel et all 



1991 . 



Bell fc de Jong! (2001) discuss mass-to-light ratios in 
the optical and near-IR passbands. Assuming an universal 
spiral galaxy initial mass function (IMF), they find vari- 
ations of the mass-to-light ratio of up to 7 in the optical, 
3 in the i3— Band and 2 in the if —Band. To convert from 
luminosities to solar masses, we us e a if —Band mass-to - 
light ratio Mk/Lk = 0.5 Mg/ Lg ljBell fc de Jonej|200ll ). 
The stellar surface density is then determined from the if - 
Band luminosities via E* = 4.27 10"^ x Lk Mq pc^^ for 
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a distance of 8.4 Mpc to M51 and an angular resolution 
of 3". 

The radial average of the stellar surface density 
is shown in Fig. [71 It decreases from a central value of 
Eh,=81OM0 pc~^ to TMqpc"^ before increasing to 12 
Mqpc~^, due to the influence of the neighboring galaxy 
NGC5195. 

The fraction of stellar-to-gaseous mass density 
S*/Sgas=10.5 is high in the center, decreasing outwards to 
values below 1 at a radial distance of about 8 kpc, before 
increasing again due to NGC5195. 

5. Kinematics of the ^^CO 2-1 data 

5.1. Velocity field and streaming motions 

Figure[9ji shows the observed velocity field of the ^^CO 2~ 
1 data ranging from 370 to 550 km with the northern 
part blue-shifted towards us and the southern part of the 
galaxy red-shifted due to the inclination of the galaxy. 

To calculate the velocity field for purely circu- 
lar motions, we u s e the rotation curve derived by 
Garcia-Burillo et all (|l993bl lah and a position angle of 
®mT" = ^"^^ °- addition, we assume a constant inclina- 
tion oi i = 20° and a central ve locity of Vsys =472 km s"^ 
(|Tu11vII1974I : IShettv et all2007t l. The velocity field at each 



point in the map is then given by 

Vmod = Vsys + v{R) X cos(e - QmD X sin(i). 

The resulting map is pre- 



(5) 



with the azimuthal angle Q. 
sented in Figure[5)3. 

The two maps resemble well, i.e. the dominant compo- 
nent creating the observed velocity field is the differential 
rotation curve v{R) creating the striking X-shape of the 
iso- velocity contours. 

However, the observed velocity field shows in addition 
many complex structures and distortions not expected 
from purely rotational motions. We define the residual 
velocity as the difference of the observed velocity field 
and purely rotational motions from the rotation curve: 
Avres = fobs — 1'mod • FigurefTUl shows a map of the resid- 
ual velocities, which we will call streaming motions hence- 
forth. 

Streaming motions vary between ±45 km s~^. Strong 
residual velocities are seen in the center and near the in- 
ner spiral arms, e.g. near (0",-85") and (60" ,20"). Both 
outer arms do not show significant streaming motions. We 
speculate that their locations at radii of ^ 120—170" from 
the galaxy center, lie close to the predicted location of the 



corotation radius at 160" ( Garcia-Burillo et al. 1993al lbl) 



At this position no net streaming is expected. 

Here, we don't go into more detail, as stream- 
ing motions in M51 have p reviously been discussed by 



Garcia-Burillo et al. (11993 



bv lAalto et al 



using older 30m data and 

, using OVRO data. In addition, 

Kuno fc Nakail (|l997l ) discussed the streaming motions in 
M51 based on ^^CO 1-0 d ata observed with th e NRG 45m 
telescope. More recently, IShettv et ah! (|2007[) use BIMA 



^^CO 1-0 data of the entire galaxy at 6" resolution to 
discuss radial and tangential streaming motions in detail. 
They determined the systemic velocity, position angle, and 
inclination from a two-dimensional fit to the data. 
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Fig. 10. Map of non-circular streaming motions Awj-es, i-e. 
the deviations of the observed velocity field from circular 
motions following the rotation curve. 
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Fig. 12. Radial averages of the velocity dispersions a = 
FWHM/(2^2 ln(2)) of the molecular gas traced by CO, 
H I, and stars. 
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Fig. 9. Left: Map of the first moment Vohs = J Tvdv/ J Tdv showing the intensity weighted mean velocity of the 
^^CO 2-1 data. Data below the 5-<7-level wer e blanked. Right: Map of the calculated velocity field Vmod using the 



rotation curve of iGarcia-Burillo et al.l (|l993bl lal). The mean velocities range from 370 to 550 km s ^. Contour levels 
are in steps of 7km s~^. 



Table 3. Equivalent widths Avoq of CO and Hi in M51 5.2.1. ^"'CO 2-1 

averaged over the areas shown in Fig [2] 





AVeq(CO) 


AVcq(Hl ) 




[km s-i] 


[km s~^] 


inner spiral arms 


25.9 


21.5 


outer spiral arms 


11.1 


18.3 


interarm region 


7.3 


19.4 


ratio inner/outer spiral 


2.3 


1.2 


inner arm/interarm 


3.6 


1.1 


outer arm/interarm 


1.5 


0.9 



5.2. Velocity dispersion of the molecular and the 
atomic gas 

The velocity dispersion of the gas is a signature of the tur- 
bulent interstellar medium and systematic motions within 
the beam. The turbulent medium and the differential rota- 
tion curve hinders the medium to collapse under its grav- 
itational pull, as described by the Toomre criterion. 



FigurelTTk shows a map of equivalent ^^CO 2-1 line widths 
Aveq = J Tdti/Tpk measuring the velocity dispersion^ 

We will refer to Aucq in the following as Awobs the 
observed line width. 

The broadest widths are observed in the center, reach- 
ing up to 100 km s^^. In the disk of M51 spiral arms show 
enhanced line widths relative to the interarm regions. The 
inner spiral shows strongly enhanced widths relative to 
the inter arm molecular gas. Along the inner spiral, the 
widths decrease from about 70 to 20 km s~^. At the same 
time, the contrast between the arm and interarm disper- 
sion diminishes. Near (60" ,100"), arm and interarm molec- 
ular gas shows the same dispersion of about 15-20 km s^^. 
Overall, regions with distances greater than aboutg 100" 
from the center show significantly lower widths of 15 to 
20 km s^^. This drop of CO line widths with radius is also 
seen in a plot of radial averages (Fig. [T2|) and in the table 
of line widths averaged over the inner, the outer, and the 
interarm regions (Table[3]). The outer south-western spiral 
arm shows uniform widths of around 15 km s^^. 



For a Gaussian line profile, the equivalent width corre- 
sponds to the full width at half maximum (FWHM). 
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Fig. 11. Velocity dispersions of the molecular and atomic gas in M51. (a) map of observed widths of ^^CO 2-1 in km 
s^^. (b) Velocity dispersion of Hi, i.e. map of the square root of the second moment in km s^^. Positions where the 
Hi drops below the Scr level of 1.33 M0 pc^^ are blanked. 



Influence of systemic motions The ^^CO 2-1 line widths 
map in Figure fTTk shows no strong systematic broadening 
of the observed line widths Auobs with varying angle in 
the disk which would indicate that a major fraction of the 
observed line widths originates in the galactic rotation of 
M51 (modeled in Figure[9|D) polluting the observed line 
widths Awobs- The velocity gradient within the observing 
beam (FWHM=11") due to the velocity field will broaden 
the observed line widths and thus lead to an overestima- 
tion of the intrinsic line width Auobs- bi the following, we 
will quantitatively estimate the magnitude of this effect 
on the observed line widths at some selected positions in 
the galaxy disk where velocity gradients are suspected to 
be the highest due to projection: the inner region (0"/0") 
and along the minor axis. Elsewhere in the disk the gra- 
dient of the velocity field is significantly lower (Figure[9}3) 
and additionally the observed line widths are similar to 
the outer parts of the minor axis (FigurefTTk) . 



near the minor axis of the galaxy from the center outwards 
(Figureini). 

The relative error of Awobs is 12% at the center and 
decreases outwards along the minor axis with values of 
less than 3% at distances 90" away from the center. This 
relative error is shown together with Augys for each of the 
six studied positions in Figure[T3l We estimate the effect of 
the polluted Awobs on Qtot (next section) using Gaussian 
error propagation for a reasonable range of values. The 
absolute change of Qtot for an overestimation of the line 
widths of 10% for the gaseous components (the broadening 
of stellar dispersion due to systemic motions is negligible) 
results in an maximum absolute uncertainty of 0.16. As 
this effect is minor, we will in the following neglect the 
broadening due to rotational motions. 



5.2.2. Hi 



Sakamotol ( 1996f ) corrected for the broadening due to 



systemic motions Awgys using: Avf^^ 
Avgys is the velocity variance across the beam as caused 
by the modeled large scale velocity field. We apply this 
method to quantify the influence of the rotational mo- 
tion at the central position (0",0") and six other positions 



Although line center of CO and Hi emission usually 
coincide quite well, dispersions can be quite different. 
Figure fTTb shows the square root of the second moment 

M2^^ = ^Jt v'^ dv/ J T dv of the Hi emission as a measure 

1 /2 

of its velocity dispersion. For a Gaussian line profile, 
corresponds to the FWHM. 
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Fig. 13. The ovcrcstimation of the intrinsic line widths 
by the observed line widths near the minor axis 
A(Auint/Auobs) in blue circles and the systemic velocity 
gradient in the beam Avgys in red triangles for the ^^CO 
2-1 data. 



The Hi line widths show values of aromrd 30-50 km 
s^^ in the central part. In general, the spiral arms do 
not show up prominently and the interarm regions show 
a larger dispersion than the spiral arms, in strong con- 
trast to the corresponding CO map. For example, the in- 
ner interarm region near (0",100") shows velocity widths 
reaching values of up to 40 km which is significantly 
higher compared to the arm regions. Towards the compan- 
ion galaxy, we see enhanced widths, comparable to those 
in the central region. H i shows a strong north-south asym- 
metry with significantly lower line widths in the souther 
parts of the galaxy. This is less pronounced in the CO- 
data. Several of these observations are reflected in the plot 
of radial averages (Fig. [T2|) and in the table of line widths 
averaged over the different regions (Table[3|). 

In the plot of radially averaged dispersions fFig. fT^ . 
the dispersion of CO gas exceeds the dispersion of Hi for 
radii up to 4kpc. The decrease of the dispersion of the 
atomic component in the inner part from around 20 to 10 
km s^^ within 2 kpc is steeper compared to the decrease of 
the CO dispersion. The radial averaged atomic component 
then stays constant at around 10 km out to radii of 
12 kpc. The CO dispersion drops to 6km s^^ below the 
dispersion of H i at radii of around 6 kpc and increases to 
11 km outwards due to the increased dispersion in the 
center of the companion. 

Our findings are consistent w i th the following sc enario 
(cf. iGarcia-Burillo et all Il993bl : iRand et all Il992l ). The 
dispersion of the molecular gas is enhanced in the spiral 
arms due to frequent collisions of molecular clouds which 
may be enhanced by streaming motions driven by the den- 
sity wave. The atomic clouds in the arms appear to form 
a different population as their local velocity dispersion is 
much less compared to the molecular gas. We speculate 
that thus Hi in the arms stems mainly from photodis- 



sociated H2 in GMCs, where the dispersion has already 
decreased due to dissipation. The interpretation of VLA 
Hi and OVRO CO observations bv lRand et aD (|l992h also 
favors this dissociation scenario. 

The plot of radially averaged dispersions (Fig. [12]) also 
shows the dispersion of disk stars. This is an important 
parameter for determining the gravitational stability of 
the disk and will be used in the following section. To es- 
timate the stellar ve locit y dispersi o n in M51, we follow 
Boissier et aTl (|2003[) and iBottemal (|l993l ). They showed 
that the velocity dispersion follows an exponential fall-off 
Cstars = CTo exp(— R/2hB) depending on the scale length 
of the B-band Iir of the disk. The B-b and scale length 
is he = 2.82kpc ( Trewhella et al. 2000l ) and the central 



stellar velocity dispersion in M51 is ctq = 113 km s 



(|McElrovlll995f ). At radii upto 10 kpc, the velocity disper- 
sion of the stars is larger than that of the molecular and 
atomic gas components. 

6. Gravitational stability 

6.1. The combined Toomre parameter Q^otai 

We investigate the gravitational stability of the disk 
of M51 using the Toomre Q-parameter, taking into ac- 
count the stellar and gaseous contribution. In general, the 
Toomre stability criterion is depending on the epicyclic 
frequency k{R), the velocity dispersion a{R) of t he com- 
pone nt considered, and the surface density E ( Toomrel 
1 1964 ). neglecting e.g. the influence of magnetic fields: 



K{R)a{R) 
ttGS ' 



(6) 



The gas should coUaps if Q drops below 1, when grav- 
itation dominates over dispersion and the epiclyclic fre- 
quency. 

We treat the molecular gas, the atomic gas, and the 
stars in the disk of M51 as three independe nt isothermal 
fluids using the expansion in wavenumber bv lWang fc Silk 
(|l994l ) and neglecting higher order terms in the wavenum- 
ber i.e. in the velocity dispersions. The model assumption 
is that the stars move through the gaseous medium with- 
out interaction. For this first order analysis, we assume in 
addition that the dense molecular clouds move frictionless 
through the diffuse atomic gas. The total Toomre param- 
eter then is: 



?tot 



(7) 



For k{R), we use the rotation curve derived from CO, 
assuming that it holds to the Hi gas and for the stars 
in the disk. For the velocity dispersion (t{R), we use the 
dispersions derived above from CO, Hi, and the stars, 
respectively. 

To calculate Q*, we used the 2Mass /f-band image 
smoothed to 11", as described before. In Fig[T31 the radial 
averages of Qtot and Qgas = {QbI + Qii ^^'^ shown. 

Taking into account only the Hi and CO gas, Qgas 
drops from values of more than 10 in the center to values 
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of between 2 and 5 for radii beyond about 1 kpc. as has al- 
ready been discussed in Paper I. In contrast to Qgas, Qtot 
is smaller and variies only little between 1.5 and 3, at all 
radii, staying always close but slighty above the critical 
value for gravitational instability of 1. This finding indi- 
cates a self-regulating disk, keeping the total Toomre Q 
parameter near 1. We will enter into the discussion of our 
results in the next section. 

The importa nce of the ste l lar co ntribution has also 
been studied by iBoissier et al. who presented ra- 

dial averages of the Toomre Q-parameter for their sample 
of galaxies, comparing Qgas and Qtot- They assume a con- 
stant velocity dispersion and find that the stellar contri- 
bution lowers the Q-parameter by up to 50% towards the 



10 



critical threshold. iLerov et al.l (|2008l ) very recently also 



studied the radial averages of the Q-parameter including 
the stellar contribution in a sample of nearby galaxies from 
THINGS. They find a marginally stable Qtot in their sam- 
ple of 23 nearby dwarf and spiral galaxies in accordance 
with our M51 Qtot-results. 

6.2. The influence of streaming motions 

Below, we will show that the contribution of non-circular 
motions to the measured velocity dispersions is small. The 
influence of streaming motions on the local gravitational 
stability can therefore be ignored. 

To test the influence of non-circular streaming motions 
on the stability, we correct in the following the observed 
velocity dispersion for these residual velocities Auios, as- 
suming Gaussian line profiles: 



(Vobs 



t^modj • 



(8) 



This reduced dispersion leads to a reduction of Qh2! 
assuming that the non-circular streaming motions do not 
stabilize the gas against collapse. 

A local change of the observed line widths Awobs prop- 
agates into Qgas according to Gaussian error propagation. 
Typical values of Qhi = 1 — 5, Qhj = 1 — 10 , aco — 10km 
s~^ and a variation of 50% of the velocity dispersion, yields 
absolute variations of Qgas of upto 56% only. 

In the radial averages, the maximum relative change 
between the corrected and uncorrected Qtot is only 7%. 
In the following, we therefore neglect the streaming mo- 
tions in the Toomre analysis, and use directly the observed 
velocity dispersion for the calculation of the Q-parameter. 

6.3. Tine local distribution of Qtot 

The local values of the Toomre Q-parameter Qtot at each 
position are studied in the following. Numerical simula- 
tions show that the local Q values c an deviate signifi- 
cantly from the radially averaged values. IWada fc NormanI 
(19991) use two- dimensional hydrodynamical simulations 
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Fig. 14. Radial average of the total Toomre parameter 
Qtot including the stellar contribution. In red the radial 
average of the gas-only Toomre parameter Qgas is shown. 
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to study the gravitational stability of the central region of 
a galactic disk. While global Q-parameter values of their 



Fig. 15. Map of the total Toomre parameter Qtot includ- 
ing the stellar contribution. We blanked pixels with Egas 
below the 3(T-level. The contours show Qtot= 0.5, 1, 2, and 
5. 

model disk are significantly above the threshold for grav- 
itational collapse, locally the disk shows unstable regions 
spatially correlated with cold and dense clumps. 

Figure fTSl shows a map of the Toomre Q-parameter 
Qtot in M51. In general, spiral arms show smaller values of 
Qtot, i-e. they are more prone to collapse under the grav- 
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itational pull, than the surrounding intcrarm gas. This is 
expected as star formation occurs predominantly in the 
spiral arms. 

The inner spiral arms show Qtot values between 1 and 
2 near to the threshold for gravitational stability. The re- 
gions between the inner spiral arms show values up to a 
factor of 2 to 3 higher as they lack surface density. The 
outer south-western spiral arm is almost everywhere criti- 
cal to gravitational collapse with Qtot below 1 or around 1. 
The outer north-eastern spiral arm shows Qtot of around 1 
and below and few small regions reaching values of about 
2. 

In summary, we see three systematic trends in Qtot : • a 
slight decrease from the central parts to the outer regions. 
In the spiral arms, peak values decrease from ~ 1.5 — 3 to 
~ 0.5 — 1. Note that in Paper I no systematic change of 
the gas depletion times from the inner to the outer parts 
was found. • Secondly, at a fixed radius, Qtot approaches 
the regime of instability in the arms, while the intcrarm 
regions are stable. • Thirdly, the outer disk beyond the 
outer spiral arms the disk is stable as well. 

A similar study at higher spatial resolution of the total 
Toomre parameter combining the stellar and gaseous con- 
tribut ion Qtot has been done in the LMC by lYang et al 



(|2007r ) using the assumption of a constant velocity disper- 
sion. As the LMC has n o prom inent spiral structure, the 
bv lYang et aP (l2007l) also shows no spiral arms 
with enhanced instability. They also stress the importance 
of the stellar contribution for studying the gravitational 
stability. Correlating the locations of massive young stellar 
obj ects (YSOs) in the LMC with the gravitational stabil- 
ity, lYang et all (l2007l ) find that 85% of the YSOs He in- 
side gravitationally unstable regions. This underlines the 
importance of gravitational stability of the disk for large- 
scale star formation in the LMC. 

Both galaxies, the LMC and M51, show large regions 
near the critical value for collapse indicating possible sites 
for star formation. The observation that Qtot is near 1 in 
large areas of the disk is notable and was fo und previously 



in ra dial averages of several other galaxie s (iBoissier et £ 



20031) . This has led several authors, e.g. ICombed (l200lli 



to speculate that self-regulation drives these galaxies near 
the threshold of gravitational stability. 

We will briefly sketch the feedback-cycle in the fol- 
lowing. The birthplaces of stars are deeply embedded in 
dense cores of molecular clouds formed through gravita- 
tional contraction and collapse of molecular material after 
loosing much of their turbulent support. Both, bipolar out- 
flows and expanding H ii-regions from young stars, stellar 
winds and supernovae explosions from evolved stars, will 
introduce mechanical energy into the interstellar medium, 
locally enhancing the dispersion of the gas, increasing the 
Q-value. In addition, ionizing radiation from young stars 
will lead to photo-dissociation of the surface regions and 
the formation of atomic gas from the molecular material. 
These feedback mechanisms are supplemented by the den- 
sity wave, enhancing cloud-cloud collisions due to orbit 
crowding in spiral arms, also leading to an increase of the 



velocity dispersion of the molecular gas inside the arms. 
Star formation in the remaining molecular gas must await 
dissipation of the enhanced macroscopic turbulence to cas- 
cade down to smaller scales and dissipate. Globally, this 
self- regulation may explain why Q stays near 1. 

In the interarm regions the gas stays stable (Qtot > 1) 
as no stars are formed and the feedback mechanisms are 
not at work. The regions critical to collapse in the outer 
spiral arms coincide with the finding of increased star for- 



mat iori efficiency in these areas by iGarcia-Burillo et al 
(1993b) based on Ha data. 



Shettv fc Ostriker ( 2008f l use hydrodynamic simula- 



tions, including feedback to study star formation and the 
structure of galactic disks. They find that a simple feed- 
back mechanism with one star formation event per cloud 
can not sustain a global spiral pattern in galactic disks. 
This indicate s that the modeling of feedback needs to 
be improved. iQuillen fc Bland-Hawthornl (j2008l ) explain 
episodic star formation in galaxies by strong feedback 
mechanisms using numerical models. 

In recent 3D smoothed par ticle hydr odynamics simu- 
lations of isolated disk galaxies. iLi et al. I {2005, 200fi) con- 
firmed the major role of gravitational stability in the pro- 
cess of star formation and were able to examine a threshold 
surface density. This evidence for the Toomre formalism 
in their simulations is also dominating any magnetic ef- 
fects on star formation. Mag netic fields slow the collapse 
of the gas but do not stop it (jShettv fc Ostriker.. 2008.) . 



7. Summary 

Using our complete ^^CO 2-1 map of M51 c ombined with 



new HI VLA data from the THINGS team (Walter et al 



20051 I2OO8I ). we created a total gas density map and a 



ratio map of molecular and atomic surface density. The 
total gas density exhibits an underlying exponential disk 
similar to the one found in 85 //m dust continuum emis- 
sion bv iMeiierink et all (|2005[ ). Maps of ^^CO 2-1 veloc- 
ity field, line widths and H i line widths were presented. 
Combining these information the combined gaseous and 
stellar Toomre Q-parameter was discussed radially and lo- 
cally in the disk of M51. The main results of this analysis 
are: 

— The averaged total gas density on the spiral arms is 
at around 30 Mq pc~^ which is a factor of 3 higher 
compared to the interarm regions. 

— The ratio of molecular to atomic gas surface density is 
highest in the inner spiral arms at around 50 decreas- 
ing to 0.2 in the outer spiral arms. This ratio shows 
a power-law dependence on the hydrostatic pressure 
-Phydro with an index of 0.87, similar to that fou nd in 
other nearby galaxies ( Blitz fc Rosolowskv 2006[ ). 

— We fit an underlying exponential disk to the total gas 
density data with a scale length of 7.6 kpc contain- 
ing 55% of the total mass. This is comparable to an 
underlying exponent ial disk fitted to the 8 50/.tm dust 
continuum map by ( Meiierink et al.l |2005( ) . The gas- 



14 



Hitschfeld et al.: Gravitational stability of M51 



dust ratio in the exponential disk stays nearly constant 
with galaetocentric radius at around 25. The radially 
averaged map of the gas-to-dust ratio in the residual 
data shows a variation between 80 and 200. 

— The ^^CO 2-1 velocity field shows distinct differences 
with respect to a purely rotational velocity field. We 
discuss a map of these residual or streaming motions. 
The observed width of the ^^CO 2-1 lines decrease 
from 100 km s~^ in the central part to less than 15 km 
s~^ in the outer regions. The widths of the Hi spectra 
show similar values, but a significantly different distri- 
bution. While the CO widths show higher dispersions 
in the spiral amrs, the H i shows the maximum line 
widths in the intcrarm regions and lower values on the 
spiral arms. 

— We estimate the gravitational stability using a com- 
bined stellar and gaseous Toomre parameter Qtot ■ The 
impact of the stellar component on the stability is sig- 
nificant reducint the Q parameter. It decreases by up 
to 70% towards the threshold for stability Qtotai = 1- 
Qtotai varies between 5 in the interarm regions de- 
creasing to values below 1 on the spiral arms. In gen- 
eral, the spiral arm regions are closer to the threshold 
for gravitational instability. The obtained values close 
to the threshold of 1 indicate self-regulation in the 
disk of M51. The self-regulation mech anism was dis- 
cussed on theoretical gr ounds by (e.g. Combed 2001 



Shettv fc Ostrikeiil2008[ ). Similar studies based on ex- 



tragalactic o bservations in M8 3 (jLundgren et al ] l2004h 
or the LMC ( Yang et al.ll2007l ) obtained also Q-values 
close to the threshold for stability. 
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